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CT  (computed  tomography)  reconstructions  of  fuel  cell  components  of  a  yet  unrivaled  spatial  resolution 
and  quality  are  presented.  This  is  achieved  by  application  of  the  novel  DIRECTT  (Direct  Iterative  Recon¬ 
struction  of  Computed  Tomography  Trajectories)  algorithm.  We  focus  on  two  different  key  issues  which 
essentially  rule  the  fuel  cell’s  durability  on  different  length  scales  and  physical  interactions.  On  the  reso¬ 
lution  scale  of  some  100  |mm  agglomerations  of  condensed  water  in  flow-field  channels  are  detected  by 
means  of  quasi -in  situ  neutron  CT  (after  operation).  Five  orders  of  magnitude  below  nanometer  sized  Ru 
catalyst  particles  on  carbon  black  support  are  visualized  by  electron  tomography.  Both  types  of  exper¬ 
iments  are  especially  adapted  to  the  type  of  material  involved  but  they  are  accompanied  by  severe 
deviations  from  ideal  CT  measuring  conditions,  as  well.  In  order  to  overcome  the  tremendous  reconstruc¬ 
tion  artifacts  of  standard  algorithms,  we  employ  DIRECTT  which  is  described  in  detail.  Comparisons  of 
DIRECTT  reconstructions  to  the  conventional  filtered  back  projection,  prove  the  significant  improvements 
in  both  experimental  methods. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

For  a  more  detailed  understanding  of  polymer  electrolyte  mem¬ 
brane  fuel  cells  (PEMFC)  very  different  information  about  their 
functional  components  have  to  be  obtained.  This  refers  to  struc¬ 
tural  details  and  ageing  behavior  which  have  impact  on  their 
operation  and  durability  [1].  The  most  direct  way  is  visualization 
in  3  dimensions.  Computed  tomography  (CT)  is  a  suitable  tool 
for  non-destructive  investigations  of  the  spatially  resolved  “den¬ 
sity”  information  from  micrometers  down  to  the  nanometer  scale. 
Different  issues  of  interest  require  sensing  tools  adapted  to  the 
respective  spatial  resolution  and  contrast  scale  to  be  resolved. 

The  fuel  cell’s  water  management  is  of  special  interest  since  it 
determines  the  gas  flux  (reactants)  in  the  flow-field  channels,  the 
degree  of  humidity  in  the  gas  diffusion  layer  (GDL)  and  the  pro¬ 
ton  exchange  membrane,  as  well  as  critical  points  of  corrosion  [2]. 
Neutrons  are  selected  for  imaging  water  agglomerates  since  they 
provide  a  rather  high  absorption  contrast  in  relation  to  the  other 
materials  [3].  Strictly  speaking,  the  contrast  does  not  arise  from  real 
absorption  but  from  scattering  effects  of  the  hydrogen  core  (pro¬ 
ton).  Here,  we  characterize  the  feasibility  of  neutron  CT  in  order  to 
distinguish  water  droplets  at  anode  and  cathode. 
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The  basic  fuel  cell  function  of  power  generation,  i.e.  the  electro¬ 
chemical  reaction  of  oxygen  and  hydrogen,  is  catalyzed  by  very 
small  metallic  particles.  Here,  electron  tomography  is  utilized 
to  image  nanometer  sized  Ru  catalyst  particles  with  the  high¬ 
est  precision  possible  in  order  to  estimate  their  free  (reactive) 
surface  which  rules  the  actual  reaction  rate.  The  outstanding 
feature  of  electron  tomography  technique  is  its  unique  spatial 
resolution  in  the  nanometer  range.  Its  first  and  most  popular 
applications  originate  from  molecular  biology  where  functional 
groups  of  biological  cell  structures  are  visualized.  However,  elec¬ 
tron  tomography  represents  a  major  challenge  to  the  robustness 
of  reconstruction  algorithms  (away  from  standard  routines)  since 
the  measured  raw  data  is  a  unique  combination  of  restrictions  (see 
experimental)  which  have  led  to  unavoidable  strong  artifacts  in 
previous  attempts  when  standard  reconstruction  algorithms  were 
applied. 

Several  CT  studies  regarding  specific  aspects  of  fuel  cells  are 
known  from  the  literature.  High  resolution  X-ray  CT  using  advanced 
X-ray  optics  was  applied  to  image  the  three-dimensional  (3D) 
microstructure  of  solid  oxides  in  SOFC  [4-6].  X-ray  microtomog¬ 
raphy  was  used  for  imaging  the  gasdiffusion  layer’s  fibre  structure 
[7,8].  Neutron  tomography  studies  aimed  at  elucidating  the  water 
management  in  complete  PEMFC  [9,10].  A  very  recent  paper 
employs  X-ray  (synchrotron)  tomography  to  visualize  the  water 
distribution  and  formation  in  the  GDL  and  flow  field  channels  of  a 
specially  designed  PEMFC  [11]. 
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Fig.  1.  The  DIRECTT  reconstruction  principle  by  matrix  plots.  The  (measured)  sinogram  (A)  is  reconstructed  (B),  a  percentage  of  B  is  selected  (C)  and  put  into  a  partial 
reconstruction  sum  (D)  which  is  projected  again  (E).  The  difference  A-E  is  the  input  for  the  next  iteration  cycle  until  the  projected  sinogram  (EN)  converges  towards  the 
initial  data  (AO). 


2.  Reconstruction  algorithm  DIRECTT 

The  DIRECTT  [12]  algorithm  represents  a  promising  alterna¬ 
tive  to  conventional  reconstruction  algorithms  such  as  FBP  (filtered 
back  projection)  or  ART  (algebraic  reconstruction  technique)  [13]. 
Several  successful  applications  and  model  calculations  have  been 
reported  elsewhere  [12,14].  At  the  current  state  of  computation 
there  are  numerous  refined  program  versions  arranged  for  the 
requirements  of  very  specific  tasks  such  as  bad  raw  data  (regard¬ 
ing  geometry,  intensity,  background)  or  severely  limited  datasets 
(deviating  from  ideal  CT  conditions).  Here  we  explain  the  basic 
ideas  and  steps  of  the  algorithm. 

The  iterative  philosophy  of  DIRECTT  is  schematically  displayed 
in  Fig.  1.  The  algorithm  is  applicable  to  parallel  as  well  as  fan  or 
cone  beam  geometry  of  projections.  Here  we  study  parallel  beam 
projections.  Fig.  1  is  consists  of  matrix  plots,  and  in  the  following 
we  refer  to  the  matrix  entries  (letter  number)  to  describe  the  single 
steps  of  reconstruction. 

We  start  with  a  so-called  density  sinogram  (AO)  which  is  the 
usual  experimental  input  after  intensity  normalization  (actually 
intensity  is  the  measured  quantity)  and  logarithmic  conversion 
according  to  Lambert-Beer’s  law  of  attenuation.  In  Fig.  1  the  hor¬ 
izontal  axis  of  the  sinograms  represents  a  detector  line  while  the 
vertical  axis  corresponds  to  the  projection  angle.  Here  we  use  com¬ 
puted  projection  data  of  the  related  project  logo.  This  corresponds 
to  the  Radon  transform  [15]  which  is  achieved  by  integration  of 
attenuation  coefficients  along  the  beam  direction  under  different 
angles  of  sample  rotation.  The  given  density  sinogram  can  option¬ 
ally  be  filtered  along  the  detector  direction.  This  is  helpful  to  avoid 
artifact  formation  equivalent  to  the  effects  of  an  unfiltered  back 
projection. 

Demanding  that  each  element  of  the  reconstruction  array  corre¬ 
sponds  to  exactly  one  sinusoidal  trajectory  of  the  sinogram  integral, 
the  values  for  all  possible  trajectories  are  computed  and  put  into 
the  reconstruction  array  (BO).  The  DIRECTT  algorithm  selects  pix¬ 
els,  i.e.  area  elements  corresponding  to  trajectories  of  dominant 
weight  ranking  within  a  predefined  top  percentage  of  all  trajec¬ 
tory  weights  for  an  update  of  the  reconstruction  (CO).  The  selected 
fraction  is  then  weighted  with  a  pre-defined  factor  and  added  to 
the  partial  reconstruction  sum  (DO)  obtained  in  the  previous  step 
(which  is  zero  in  the  first  cycle). 


An  intriguing  feature  of  DIRECTT  refers  to  the  possibility  to  adapt 
the  filter  function  for  trajectories  evaluations  of  specific  interest 
such  as  mass  or  contrast.  The  weights  of  all  possible  trajectories 
(corresponding  to  integer  reconstruction  positions)  are  computed 
by  averaging  along  the  respective  traces  within  the  optionally  fil¬ 
tered  sinogram.  This  involves  interpolation  between  neighboring 
sinogram  pixels  corresponding  to  different  detector  elements.  This 
is  a  substantial  difference  to  FBP,  where  interpolation  is  performed 
in  Fourier  space  (fast  FBP)  or  in  the  backprojection  step  after  the 
sinogram  has  been  (inversely)  Fourier  transformed  and  filtered 
(real  space  FBP). 

The  projection  EO  (Radon  transform)  of  the  reconstruction  array 
(i.e.  a  computed  sinogram)  is  then  subtracted  from  the  original 
data  set  (AO).  The  obtained  residual  sinogram  (Al  =  AO  -  EO),  con¬ 
taining  trajectories  of  the  remaining  elements  is  subject  to  the 
same  procedure  in  the  subsequent  iteration  steps.  Now  the  par¬ 
tial  reconstruction  sum  D1  is  calculated  from  DO  +  Cl ,  or  in  general 
Dx  =  Dx  -  1  +  Cx,  where  x  is  the  number  of  iterations.  This  is  contin¬ 
ued  until  a  pre-selected  criterion  of  convergence  is  reached.  Under 
ideal  conditions  the  residual  sinogram  converges  to  zero  which 
means  that  the  measured  (or  simulated  sinogram  AO)  is  identical 
to  the  computed  projection  EN  of  the  last  partial  sum  DN.  There¬ 
fore  the  reconstructed  structure  corresponds  to  the  spatial  density 
distribution  of  the  projected  sample. 

The  overall  procedure  corresponds  to  an  iterative  Radon  trans¬ 
form.  In  contrast  to  FBP  there  is  no  integral  computation  along 
the  line  detector  (including  limited  sampling  due  to  its  element 
size)  but  an  optional  oversampling  along  the  numerous  projec¬ 
tion  angles.  One  of  DIRECTT’s  unique  characteristics  is  its  very 
precise  projection  of  reconstruction  elements  taking  into  account 
their  actual  size  and  shape  which  is  essential  for  enhanced  spa¬ 
tial  resolution.  That  is,  reconstruction  pixels  are  considered  as  a  set 
of  densely  packed  elements  instead  of  being  (circularly  smeared) 
point  functions  only.  Hence,  all  previously  described  calculations 
are  performed  based  on  squared  area  elements  in  a  Cartesian 
matrix. 

DIRECTT  is  of  particular  interest  when  the  focus  is  on  recon¬ 
struction  of  finely  structured  details  or  on  precise  location  of 
reconstructed  elements  rather  than  on  computing  time.  In  con¬ 
trast  to  FBP,  the  DIRECTT  procedure  does  not  treat  each  (detector) 
projection  individually  by  filtering  of  any  kind.  Instead  the  entire 
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Fig.  2.  Selected  neutron  CT  density  sinogram  of  threefold  fuel  cell  stack  (top);  bottom:  graphs  of  horizontal  cumulative  cross  section  (bold  line)  and  linear  cross  sections  at 
the  detector  edges  (thin  lines). 


trajectory  of  a  reconstruction  element  is  considered  over  all  pro¬ 
jection  angles.  In  contrast  to  ART,  DIRECTT  does  not  modify  all 
reconstruction  elements  simultaneously. 


3.  Experimental 

3.1.  Neutron  CT 

Neutron  tomography  experiments  have  been  performed  at 
beamline  CONRAD  at  Helmholtz  Centre  Berlin  (HZB,  research  reac¬ 
tor  BER II).  A  complete  threefold  fuel  cell  stack  (inch  metal  housing 
and  supply  lines)  was  run  on  site  and  investigated  in  two  states: 
immediately  after  operation  and  30  min  later,  i.e.  in  a  quasi¬ 
stationary  state  (as  described  in  [10]).  The  cell  stack  was  subject  to  a 
full  rotation  through  a  360  deg  sector.  1000  single  projections  were 
recorded  using  a  LiF  flat  panel  detector  (2048  x  2048  pixel,  100  [xm 
pixel  size,  dynamic  range  1 4  bit).  Additionally,  1 0  flat  field  and  dark 
field  images  were  collected  for  the  purpose  of  data  correction. 

Fig.  2,  top  displays  an  arbitrarily  chosen  sinogram.  A  sinogram 
is  obtained  from  a  fixed  detector  line  (perpendicular  to  the  rota¬ 
tion  axis)  at  varying  angles  from  a  pile  of  projection  images.  Fig.  2, 
bottom  shows  the  according  cumulative  cross  section  along  the 
projection  angle  (fat  black  line  in  the  graph  below).  Two  orthog¬ 
onal  single  projections  at  Odeg  (or  180  deg;  the  cell  is  irradiated 
perpendicular  to  the  single  cells)  and  90  deg  (irradiation  parallel 
to  the  cells)  are  inserted.  Ideal  CT  conditions  require  that  the  pro¬ 
jected  mass  is  a  conserved  quantity  under  all  projection  angles. 
As  to  be  seen  from  the  right  insert  the  fuel  cell  is  only  partially 
projected  under  0  and  180  deg  (and  neighboring  sections).  For  fur¬ 
ther  rotation  the  projected  mass  increases  sinusoidal.  When  the 
entire  sample  is  projected  on  the  detector  (see  the  left  insert)  a  con¬ 
stant  level  of  the  integrated  density  (mass)  should  be  expected.  But 
because  of  very  high  absorption  the  disadvantageous  aspect  ratio 
causes  nonlinear  phase  dependent  attenuation.  Therefore  only  the 
angular  section  outside  the  shaded  region  is  used  for  the  recon¬ 
structions. 

An  even  more  severe  problem  is  implied  by  the  large  neutron 
aperture  which  amounts  to  some  centimeters  and  results  to  inher¬ 
ently  un-sharp  projection  images.  This  effect  will  be  designated  as 


‘smearing’  and  will  turn  out  to  be  a  function  of  the  sample-detector- 
distance:  the  larger  distance  between  sample  and  detector  the 
more  blurred  is  the  obtained  image.  Basic  calculations  how  to 
properly  take  into  account  a  position  dependent  blurring  to  the 
reconstruction  algorithm  have  been  given  [16].  A  first  approach 
aims  at  computing  the  smearing  in  the  projection  part  of  the  iter¬ 
ative  DIRECTT  algorithm  and  neglecting  it  in  the  reconstruction 
part. 

Tracing  single  sinusoidal  trajectories  (of  marked  objects)  reveals 
a  rotation  axis  which  is  inclined  by  approximately  1 .3  deg  with 
respect  to  the  Cartesian  detector  coordinates.  Thus,  planar  sino¬ 
grams  are  oblique  planes  in  projection  space  (i.e.  the  3D  image 
stack  of  projections).  The  planar  information  is  separated  by  simple 
rotation  of  the  single  projections.  Additionally,  there  are  indications 
for  axis  misalignment  in  the  direction  of  the  (nearly  parallel)  beam 
which  yields  elliptical  trajectories  in  the  projection  space.  This  mis¬ 
match  cannot  be  corrected  by  pure  data  pre-processing  but  is  a 
challenge  for  experimental  alignment. 


3.2.  Electron  tomography 

Electron  microscopic  projections  of  catalysts  on  carbon  black 
support  have  been  performed  ex  situ  in  a  Zeiss  Libra®  Microscope. 
The  bright  field  images  (in  zero  loss  modus)  are  recorded  in  par¬ 
allel  beam  geometry  using  200  keV  electrons.  The  actual  detector 
resolution  used  here  is  0.26  nm.  Since  none  of  the  state-of-the-art 
manipulators  can  warrant  sufficient  precision  in  the  nanometer 
range  the  technique  suffers  from  an  inherently  unstable  rotation 
axis.  The  measurements  are  performed  on  samples  positioned  on 
copper  grids  with  a  thin  carbon  foil  containing  gold  markers.  In 
order  to  have  available  the  possibility  of  re-orientation  of  fluctuat¬ 
ing  single  projections  spherical  gold  markers  are  introduced  to  the 
measured  volume. 

Due  to  the  sample  holder’s  size  (relative  to  the  magnetic  pole 
shoe  distance)  a  full  rotation  of  the  sample  is  impossible,  resulting 
in  an  integral  Limited  View  problem.  The  measured  angular  section 
typically  amounts  to  140  deg.  Because  of  the  high  electron  beam 
brilliance  local  material  degradations  cannot  be  excluded.  Thus,  the 
recording  time  per  projection  (about  90s)  has  to  be  kept  as  short 
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Fig.  3.  Neutron  CT  reconstructions  of  transversal  section  of  the  threefold  fuel  cell;  left  and  right:  separate  reconstructions  of  the  two  180  deg  sectors  of  the  sinogram  applying 
a  distance  related  linear  ramp  factor  in  direction  of  each  projection;  middle:  combination  of  both  sectors. 


as  possible.  The  total  tomographic  measurement  comprises  merely 
143  projections,  i.e.  with  respect  to  the  detector  size  ( 1 550  pixel  per 
line)  the  reconstruction  suffers  from  a  tremendously  underdeter¬ 
mined  Few  Angles  problem. 

The  sample  (incl.  Au  markers)  is  much  larger  than  the  pro¬ 
jected  cross  section,  i.e.  volume  elements  near  the  edge  of  the 
reconstructed  region  are  not  projected  at  all  rotation  angles  and, 
worse,  the  integral  mass  is  not  conserved.  This  corresponds  to  the 
well-known  region  of  interest  (ROI)  problem  which  hampers  the 
empirical  determination  of  the  exact  axis  position  since  tracking 
the  centre  of  gravity  as  a  function  of  phase  is  pointless.  However, 
knowing  the  axis  position  with  the  highest  precision  possible  is  of 
the  utmost  significance  for  the  achievable  reconstruction  quality. 

High  density  volume  elements  (Au  markers,  catalyst  parti¬ 
cles)  are  nearly  opaque.  The  as-measured  intensity  is  out  of  the 
valid  dynamic  range.  As  a  result,  the  reconstructed  array  exhibits 
intense  streak  artifacts,  since  the  actually  penetrated  lengths  of  the 
absorbing  objects  cannot  be  converted  to  meaningful  values  of  the 
respective  attenuation  coefficients.  Thus  for  the  concealed  partial 
volume  a  partial  phase  dependent  Limited  View  problem  results. 

From  the  projected  sector  and  the  number  of  projections  a  mean 
angular  increment  of  1  deg  is  derived.  However  in  the  electron 
microscope  the  sample  rotation  can  only  be  pre-selected  at  an  accu¬ 
racy  of  about  0.2  deg.  The  resulting  sequence  of  non-equidistant 
projection  angles  is  created  by  alignment  and  rectification  of  the 
single  projections.  Therefore  the  application  of  an  arbitrary  list  of 
projection  angles  has  been  inserted  to  the  DIRECTT  reconstruction 
program  in  order  be  more  precise. 


4.  Results 

4.1.  Water  agglomerations  in  flow-fields 

Part  of  the  experimental  drawbacks  of  neutron  radiographs 
which  refer  to  the  phase  dependent  integral  attenuation  and  the 
suppression  of  artifacts  due  to  incomplete  trajectories  of  projected 
elements  (ROI)  have  not  been  overcome  significantly  by  iterative 
DIRECTT  reconstructions  and  did  not  justify  the  required  additional 
computing  time. 

However,  computation  of  the  position  dependent  smearing  of 
projections  (due  to  the  large  aperture)  yields  significantly  improve¬ 
ments.  It  is  applied  to  the  projection  part  of  iterations  exclusively 
as  desmearing  by  deconvolution  of  the  reconstructions  would  not 
provide  a  position  dependent  procedure  [7].  Analogously  to  the 
magnification  function  in  cone  beam  geometry,  the  smear  width 
increases  as  a  function  of  the  sample’s  detector  distance.  In  order 
to  get  a  quantitative  measure  of  the  smearing  function  calibration 
measurements  using  a  thin  Gd  foil  have  been  performed.  In  the 
result  we  obtain  a  smear  width  calibration  as  a  second  order  poly¬ 
nomial  dependency  on  the  detector  distance  including  a  4  pixel 
offset  (according  to  inherent  detector  smearing).  On  the  other  hand 
the  smear  width  is  no  function  of  the  lateral  position. 


Although  180  deg  parallel  beam  data  sets  generally  contain  all 
information  necessary  for  reconstruction,  the  full  rotation  data  is 
required  for  the  position  dependent  iterative  smearing  procedure 
in  order  to  perform  de-smearing  after  all,  thus  improving  the  recon¬ 
struction  quality.  Fig.  3  illustrates  the  successful  application  of  the 
position  dependent  de-smearing  by  simply  weighting  the  ampli¬ 
tudes  as  ramp  functions  along  the  direction  of  projection  at  the 
example  of  180  deg  sections  of  the  sinogram  (as  sketched  in  the 
inserts),  in  which  due  regions  (close  to  the  detector)  are  favored  by 
‘sharp’  projection,  respectively. 

Comparing  reconstructions  of  complete  data  sets  an  unam¬ 
biguous  advantage  of  the  DIRECTT  algorithm  using  the  position 
dependent  smearing  is  obvious.  This  is  illustrated  at  the  example 
of  the  middle  cathode  in  Fig.  4  (perpendicular  to  the  reconstructed 
slices),  where  the  differences  are  very  clearly.  The  water  droplets 
are  localized  very  well  in  the  DIRECTT  reconstruction  but  poorly  in 
filtered  back-projection. 

As  to  be  expected  the  more  significant  water  indications  (judg¬ 
ing  from  a  contrast  point  of  view)  are  found  at  the  cathode  sites 
(where  the  protons  migrate,  i.e.  where  the  water  is  formed).  In  the 
FBP  the  cathode  structures  are  partially  superimposed  by  neighbor¬ 
ing  anodes  (Fig.  4,  crossing  channels).  Obviously  the  determination 
of  the  condensed  water  concentrations  at  the  different  electrodes 
has  been  improved  by  the  DIRECTT  reconstructions. 

4.2.  Nanometer  resolved  catalyst  particles 

As  outlined  in  Section  3.2  the  non-equidistant  angular  sequence 
of  electron  tomographic  projections  has  been  applied.  Two  indica¬ 
tions  confirm  the  advantages  over  approximated  increments:  the 
projected  trajectories  of  the  gold  markers  are  much  more  sinusoidal 
and  the  reconstruction  quality  increases.  This  includes  both  local¬ 
ized  details  as  well  as  integral  criteria  such  as  weight  and  variance 
of  the  residual  sinogram. 

The  reconstruction  of  the  rather  limited  raw  data  set  by  stan¬ 
dard  FBP  results  is  very  disappointing  images.  Careful  geometry 
corrections  are  required  (e.g.  precise  axis  position)  and  an  inten¬ 
sity  threshold  (strong  absorbers)  in  each  case.  Neglecting  such  data 
pre-processing  results  in  reconstructions  dominated  by  streak  arti¬ 
facts  which  conceal  the  structure  of  interest  (Fig.  5,  FBP).  Moreover, 
these  adjustments  are  a  pre-requisite  to  enable  the  iteration. 

As  a  result  of  the  simultaneous  occurrence  of  strong  (gold) 
absorbers  and  the  ROI  problem  large  masses  (high  densities)  ‘from 
outside’  (not  reconstructed)  density  produce  strong  artifacts  inside 
the  ROI.  Within  DIRECTT  the  ‘disturbing’  mass  can  be  suppressed 
by  extending  the  reconstruction  array  and  thus  localizing  it  outside 
the  ROI.  This  is  only  possible  as  of  another  advantage  of  DIRECTT: 
the  geometrically  true  evaluation  of  incomplete  trajectories  of  mass 
elements  from  outside  the  ROI. 

The  quality  differences  for  imaging  the  few  nanometer  catalyst 
particles  by  different  reconstruction  algorithms  are  demonstrated 
in  Fig.  5.  At  the  example  of  an  arbitrarily  chosen  detail  the  FBP  of 
the  raw  data  provides  mainly  artifacts.  SIRT,  the  commonly  used 
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Fig.  4.  Comparison  of  FBP  filtered  back-projection  (left)  and  DIRECTT  (right)  reconstructions  at  the  middle  cathode  plane:  different  contrasts  and  improved  resolutions 
(perceptibility)  of  water  droplets. 


Fig.  5.  Quality  differences  of  different  reconstruction  algorithms  as  revealed  by  imaging  properties  of  catalyst  particles  on  selected  carbon  black  support  (reconstruction 
detail). 


algorithm  for  electron  tomography,  an  iterative  ART  procedure  per¬ 
forms  better  when  appropriate  data  preprocessing  is  applied.  The 
higher  quality  of  the  DIRECTT  reconstruction  is  obvious  (Fig.  5, 
right),  although  the  required  computing  time  is  two  orders  of  mag¬ 
nitude  shorter. 


Fig.  6.  3D  visualization  of  Ru  catalyst  nanoparticles  and  zoomed  detail  at  about  1  nm 
spatial  resolution:  dominant  3-8  nm  Au  markers  on  cloudy  carbon  black  supports. 


Beyond  the  2D  details  of  Fig.  5  the  3D  visualization  of  Fig.  6  gives 
an  impression  of  the  relative  sizes  and  the  catalyst  particle  arrange¬ 
ment  on  the  membrane  and  the  carbon  supports.  Only  such  precise 
results  permits  advanced  data  analysis  which  aim  at  the  estimation 
of  the  catalysts  free  surface.  This  includes  the  individual  ellipsoidal 
particle  deformations  and  the  determination  of  the  catalyst  fraction 
embedded  in  carbon. 

5.  Conclusion 

The  unique  tomographic  algorithm  DIRECTT  has  been  applied  to 
solve  difficult  reconstruction  tasks  regarding  fuel  cell  components. 
It  copes  with  aperture  smearing  effects  of  neutron  CT  and  resolves 
separately  the  water  agglomerates  at  fuel  cell  cathode  and  anode 
sites.  In  case  of  the  simultaneous  occurrence  of  numerous  restric¬ 
tions  inherent  to  electron  tomography  it  provides  advantages  over 
the  other  relevant  algorithms.  Beyond  the  application  to  fuel  cell 
components  the  reconstruction  of  any  restricted  tomographic  data 
set  by  DIRECTT  is  expected  to  yield  improved  results  compared 
to  the  state  of  the  art  provided  that  proper  data  preprocessing  is 
applicable  or  precise  measurements  are  given. 
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